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Abstract 



In a previous paper we described a quantum algorithm to prepare an arbitrary 
state of a quantum register with arbitrary fidehty. Here we present an alternative 
algorithm which uses a small number of quantum oracles encoding the most sig- 
■ nificant bits of the absolute value of the complex amplitudes, and a small number 

. of oracles encoding the most significant bits of the phases. The algorithm given 

Oh! here is considerably simpler than the one described in 01 , on the assumption that a 

sufficient amount of knowledge about the distribution of the absolute values of the 
^ ■ complex amplitudes is available. 

^3 ■ 
^■ 

1 Overview 

><■ 

d ' The first step of many quantum computer algorithms is the preparation of a quantum 

register in a simple initial state, e.g., the equal superposition of all computational basis 
states. Some applications of quantum computers, such as the simulation of a physical 
system jSJ El E] , require the initial preparation of more general states. Here we consider 
the state preparation problem in the case that the Hilbert-space dimension of the quantum 
register is so large that listing the complex coefficients of the state is impractical. 

In a previous publication we have shown how to use elements of Grover's algorithm 
[S] to prepare a register of logj qubits, with arbitrary fidelity, in an approximation to 
the state 

N-l 



=\\/p{x)e'^"f"-^^\x) (1) 



for any probabilities p{x) and phases We assume that is an integer power of 2. 

Here and throughout the paper, |0), . . . denote computational basis states. 

We will now describe an alternative algorithm to achieve the same goal, i.e., to prepare 
the quantum register in a state |\E') such that the fidelity, K^'l^')!, is close to 1. In the 
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description of the algorithm below, we introduce three positive integer parameters, T, T', 
and a. We will indicate how to choose these parameters, and derive a lower bound on 
the fidelity |(\E'|\E')| in terms of them. We will also show how the required computational 
resources scale with these parameters. This will put us in a position to compare the 
two versions of the algorithm. See Ref. [T] for a comparison with the state-preparation 
algorithms by Kaye and Mosca [H] and Grover and Rudolph [?]. 

The functions p{x) and (j){x) are assumed to be given in the form of classical algorithms. 
The function p{x) is used to construct a set of quantum oracles as follows. Let T be a 
positive integer. For = 1, . . . , T, we define 

Ok{x) = Ck{x) , (2) 

where Ck{x) G {0, 1} are the coefficients in the binary expansion 

oo 

^/^iN^=J2ck{x)2-' , (3) 

k=l 

and where is a positive real number, r] < 1, such that 

p{x) < — ^ for all X . (4) 

We extend this definition beyond the domain of the function p by setting Ok{x) = for 
X > N. The quantum oracles are unitary operators defined by 

6,|x) = (-lf'=(^)|x), (5) 

which can be realized as quantum gate sequences using the classical algorithm to compute 
the probabilities p{x). For each oracle, we define the number of solutions 



Nk = J20k{x). (6) 



Now let T' be a positive integer, and let c'j^{x) G {0, 1} be the coefficients in the binary 
expansion 



(p[X) 

k=l 

For A; = 1, . . . , T', we define unitary operations, Ui, . . . , Ut', on our quantum register by 

The operators Uk are conditional phase shifts that can be realized as quantum gate se- 
quences using the classical algorithm to compute the phases jSj. 

The algorithm can now be described as follows. Choose a suitable (small) number, a, of 
auxiliary qubits, and define L = log2 N + a. Prepare a register of L qubits in the state 

|^°) = (2'^Ar)-i/2 J2 \x) . (9) 

x=0 
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For = 1, . . . , T, define the Grover operator 

G(0,,4)= ((2|vI>°)(vI/V/)0.)*' , (10) 

where I is the L-qubit identity operator, and where the integer "times" tk are defined 
in Eq. below. Apply the Grover operators successively to the register to create the 
state 

|*^) = G(OT,tT)---G(Oi,ti)|^°). (11) 

Now measure the a auxiliary qubits in the computational basis. If one of the outcomes 
is 1 (the probability for this will be shown to be small), this stage of the algorithm has 
failed, and one has to start again by preparing the register in the state |^°) as in Eq. Q. 
Otherwise, i.e., if all a outcomes are 0, this stage of algorithm has succeeded, and the 
resulting state of the remaining L — a = logg N qubits, which we denote by \'^r), "will be 
a good approximation to the real-amplitude state 

N-l 

l^r) = $^VpMk)> (12) 

x=0 

obtained from our target state |\1') by setting the phases 0(x) to zero. 

The final stage of the algorithm adds phases to the state |\E'r) by applying the operators 
Ui,...,Ut', 

\^) = UiU2---Ut'\^,) . (13) 

In the next section, we analyse the dependence of the state |\E'"^) on the numbers tk, and 
thus motivate the definition (jl^ . At the end of the section, we derive upper bounds on 
the numbers tk and therefore on the required number of oracle calls. In the final section, 
we derive a lower bound on the fidelity in terms of the parameters T, T' and a. 



2 Number of oracle calls 

In the following we will use the notation 1 : ri to index an ordered sequence of n symbols, 
for example, 

qi.n = qi, . . . , Qn ■ (14) 
Using this notation, the statement gi:„ = Ci-n means that qj = Cj for any j = 1, . . . ,n. 
We define a set of refined oracles, 

n _ / 1 if = ci;„(x) , , . 

"^^^"^^ \ otherwise, ^ ^ 

which can be expressed in terms of the oracles Ok as follows. 

n 

0,,Jx) = l[\Ok{x)-l + qk\. (16) 

k=l 
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Let fiq^.j. be the set of values that are accepted by the oracle 0^^.^., i.e. 

^,.. = {^--0,Ux) = l}. (17) 
Furthermore, denote by N^^^ the size of the set f^^^.j., 

X 

The first stage of our algorithm takes the initial state through a series of intermediate 
states, to the state \^'^). Due to the properties of the Grover operators, G{Ok,tk), 
the intermediate states are of the form 

where 

k 

<.=B' + Y.^,h,, (20) 

i=i 

where the features hj are positive numbers determined by the times tk, and where the B'^ 
are real numbers determined by the normalization conditions (vI/'=|vl/'=^ = i. 

We show next how the features hk depend on the numbers of Grover iterations tk- 
2.1 General oracle 

We will be using the following result of Biham et. al. jH]. Consider an oracle O, which 
accepts r values (out of the total of 2°iV, i.e., Yl^x=o~^ 0{x) = r). We shall call such values 
of X good, as opposed to had values of x that are rejected by the oracle. Using different 
notation for the coefficients of good and bad states, we have that after t Grover iterations 
an arbitrary quantum state 

I*"") = E ^"'1^) + E (21) 

good X bad x 

is transformed into 

\^^^) = G{o,twn = J2 9^1^) + E ^M^)- (22) 

good X bad x 

Let g™^ and 6'" be the averages of the initial amplitudes of the good and the bad states 
respectively: 



^^" = ^E^"' ^" = ^^E^"' (23) 

good X bad x 

and similarly for the final amplitudes 

^''^ = ^E^^' ^'"^ = ^^E^'^ (24) 

good X bad x 
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Let us also define 

Agf = gT - , = - b"'' . (25) 



mi 

X 



In other words, A^fjf' and A6™^ define the features of the initial amplitude functions g 
and relative to their averages and Biham et. al. have shown that the change 
of the amplitudes is essentially determined by the change of the averages: 



= /" + A(7r 
e = l^^' + i-lfMr, (26) 



where the averages g^"^ and are given as follows. Define 



uj = arccos ( 1 ) , (27) 



a = + |^i°i|2 , (2^ 



arctan(H— / ) . (29) 



The averages are given by 



g^"" = \l a sm{ujt + (j)) 



a cos{ujt + (j)) . (30) 



These formulas allow us to calculate the number of Grover iterations t from the ratios 
^'"^76'"^' and g^'^/W'' as follows. From Eqs. ^ we have 



|- = ^^^tan(^t + 0), (31) 



which, together with Eq. (j29j) . gives 



u!t = arctan ( ^ / — arctan ( \ / . (32) 

2.2 Formulas for tk 

Consider the state i.e. the state that results after building the first k features using 
the oracles Oi, . . . , Ok- Let g^^+i^ ^fc+i t>e the average amplitudes of the "good" and "bad" 
states within l^^'^) with respect to the oracle Ok+i- By direct calculation we have 

-mi _ ^gi:fc Ql-.k gi:fc-l- 7im _ ^qi:k Ql-.k Ql:k^ ( 'i.'X\ 

^Qi-.k gi:fcl ^Qv.k Qi-.kO 

and therefore 
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Similarly, in the case of the final averages g^'^^ and b^\i we obtain 

Below we need expressions for the ratios ^^V^™' (7f°/&fc°, which follow by substituting 
k for k + 1. The number of Grover iterations, tk, required for converting the state 
into |\E''^) can then be obtained from Eq. (jS21), 



where 



Uktk = arctan I fl-i/ ^ ^ | — arctan ( f^-t/ ^ ^ 1 , (36) 



ojfc = arccos ( 1 - ^ ) . (37) 



Of course these formulas for the integer times tk are useless by themselves, because they 
depend on the coefficients A^_^ ^ , which are defined in terms of the unknown features h 



k 



[see Eq. (r2(Jj) ]. The following argument leads to an explicit formula for the tk- 
By construction of the sets ^^g^.^., the sums Yl^=i Q'i^^V fc-bit approximations to 
the target amplitudes \/ p{x) for all x G f^g^.^- We thus aim for the features hj to be as 
close as possible to the values / yJrjN . This motivates the following choice for the tk- 
Instead of the amplitudes (1201), define 

k 

where the hj have been replaced by / ^/r]N , and where the terms S'*^ are determined 
by the normalization conditions (^!''^\^)'^) = 1 for the states 

i^") = E E <J-)- (39) 

These states can be regarded as fc-bit approximations to the intermediate states l^'^). We 
thus get the following modified expressions for the average amphtudes. 



and 

The final expression for the times tk is then 



(40) 



(41) 



tk 



7^ + — h 




Nk 

„ , , arctan I — ^ \ I — — — | — arctan 

2 ujky \ h'^"" V 2-N - 




(42) 



6 



where the extra term 1/2 combined with the [. . .J operation amounts to a rounding to 
the nearest integer. 

The expressions for depend exphcitly on the numbers Ng_^^, i.e., the numbers of points 
X for which the k most significant bits of \JrjNp{x) are given by qi-k. If this global 
information about the probabihties p{x) is available, the version of our state preparation 
algorithm described here will often be simpler than the original version of the algorithm 
described in Ref . [2| . If the numbers A*"^,^ ^ are not available initially, they can be obtained 
via the quantum counting algorithm jTOj. In this case, the algorithm described here loses 
much of its appeal. The analysis of the original algorithm in Ref. includes bounds for 
the resources required for the initial quantum counting step. 



2.3 Bound on the number of oracle calls 

We have, by definition. 



2^N 2 

Since x"^ > sin^ x we obtain 



1 — coscufc = 2 sin — . (43) 



Furthermore, we have 
and hence 



^^>2\/^- (44) 



Uktk < 2tx , (45) 



tk< — <nJ--. 46 
V Nk 

The overall number of oracle calls is therefore bounded by the expression T'+Tn ^J2°'N / N^. 
A typical value for the fraction N^/N is 1/2. The worst case for the number of oracle 
calls corresponds to A^^ = 1, which is equivalent to Grover database search The ef- 
ficiency of our algorithm can be improved by ignoring very small values of A^^. Bounds 
for the corresponding fidelity reduction have been derived in Ref. [T]. An analysis of the 
asymptotic number of oracle calls in the limit of large A^ is possible, e.g., for a sequence 
of states for which the parameter t] does not depend on A^ and the ratios Nk/N tend to a 
constant Cfc as A^ — >• oo. In this case, the fidelity bound does not depend on A^. For 
the right-hand side of the bound , we have then 



2'^A^ 

VTA / ^ ^ 7rv/2VCfc as N oo , (47) 
i.e., the bound for the required number of oracle calls tends to a constant for large A^. 



3 Fidelity analysis 

In this section we derive a lower bound for the fidelity |(\E'|\E')|. We start by considering 
the fidelity between the real-amplitude target state |\l/r) defined in Eq. flT^ and the state 
|\E'^) resulting from the Grover iterations, but before the a auxiliary qubits have been 
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measured [see Eq. ffTT|) ]. It follows from the discussion at the start of Sec. |21that |\E'"^) 
can be written in the form 



(45 



x=0 



The first step is done in subsection 13. ![ where we show that the choice Eq. ()42|) for the 
integer times tk implies that the features hk-, for = 1, . . . ,T, satisfy the inequalities 



hk - 



< 



2l-a/2 



(49) 



Subsection 13.21 uses this result to derive the fidelity bound 



x=0 
N-1 



^V^(i?^ + 5^c,(x)/., 

x=0 j=l 



> 1 -3T 



2-a/2 

V 



(50) 



where we have used the fact that p{x) = for x > N, and where we have assumed that 
T is chosen to be the smallest integer for which 



2T2 



< 2" 



(51) 



One can of course use bigger values of T, but this would not improve the performance 
of the algorithm as the fidelity of the state preparation is limited by the choice of a [see 
Eq. (gni)]. 

The next step of the algorithm is the measurement of the auxiliary qubits. The probability 
of failure, pfaii, i.e. the probability of obtaining a nonzero result, is given by 



Pi,n = {2^N-N)\B 



T|2 



(52) 



where is the normalization term in the expression (j48j) for |^"^). Subsection 13 .21 derives 
the following bound on the failure probability. 



Pfail < 16T 



2-a/2 



V 



(53) 



If there is no failure, i.e., if the measurement outcome is zero, the post-measurement state 
is given by 



^ N-1 T 



Together with Eqs. (j48p and (j5Uj) it follows directly that 

1 



> 1 - 3T 



2-a/2 



(54) 



(55) 



Finally, in subsection 13.31 we combine this bound with a simple analysis of the last stage 
in which the phases are added to the real amplitudes of the state |^r)- The result is the 
following overall lower bound on the fidelity between the target state and the state 
|\E') prepared by the algorithm, 

I (^1^) \> (^1-3T (1 - 2-2^'-^) . (56) 

This bound determines the performance of the state preparation algorithm described in 
this paper. 



3.1 Upper bound on hk — 2 ^ / ^/rjN 



Consider the development of a single feature, h, in t Grover iterations based on an oracle O. 
Let r be the number of good states, or solutions, of O. It follows from Eqs. (IHUI) that h 
depends on t via 

hit) = a^f¥Wfr sin(cjt - , (57) 

where the values of a and ^ depend on the initial average amplitudes and of the 
good and the bad states with respect to the oracle O. According to Eq. (j2Hl) "we have 



a 



Mni|2 I I— ini|2 



2" - r/iV 

The average amplitude of "bad" states is bounded as 

1 



(5^ 



(59) 



V2^ 

and the maximum possible value of the average amplitude of "good" states is bounded as 

1 
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(60) 



Hence Eq. 



implies 



< 

< 
< 



1 



r/N 



2"iV r/A^(2" - r/N) 
1 1 

2^ ^ 7]N{2'' - 1) 
1 2 



4 



(61) 



Since in our algorithm the value of tk is rounded to the nearest integer, hk = h(tk) will 
rarely coincide with the target value of 2~''/ y/rjN . The mistake, however, can be bounded 
as 

hk ^ <max\h{t + l) -h{t)\, (62) 
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where the maximum is taken with respect to the quantities a, tu, ^, r and t. The param- 
eters characterizing the algorithm, a, rj and A^, are being kept constant. Using we 
obtain 



hi 



< max ■ 



rir 
2 



sin (^(tut — ^) + Lu) — sm{ujt — ^) 



< max sin(ci;/2) 

Vvr 



(63) 



where the last inequality follows from the properties of the sin function. Since < u < it 
we have 



sin(c<j/2) 



1 — cos UJ 



2«iV 



which implies the bound 



as required. 



hu - 



< 



2l~a/2 



(64) 



(65) 



3.2 Lower bound on |(^, |*^ 
Directly from the definitions we have 



x=0 
N-1 



T 



x=Q \ 7=1 * ' 1=1 ^ V / / 



x=0 \ j=l ■ ■ ] 

where we have used the fact that p{x) = for x > N. Let us define 

T 



X) 



j=T+l 



2-J 



2-3 
/r]N 



Since 



2-3 



\/ p{x) — h{x) 



and since Y1!1=q Pi^) = 1 we have 



K^rl^Ol > 1 



N~l 



x=0 



N-1 



a;=0 



Using the bound (jEHI), one can show that 



\6{x)\ < 2T 



2-a/2 



(66) 



(67) 



(6^ 



(69) 



(70) 
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and therefore, using ^/p{x) < l/^/rfN, 



>1 



N-1 



x=0 



2T- 



2-a/2 



(71) 



The function 6(2;) satisfies the bounds 

< < 2"^/v/r/iV. (72) 

In order to find the lower bound on | (\E'r|\E'-^) | from Eq. ()71|) we need to calculate \B'^\. 
This can be done by examining the normalization condition {^"^Y^"^) = 1 which reads 

2 



1. 



(73) 



x=0 



Using the definitions (jEZI), this can be rewritten as 

J2 (b^ + Vp{^ - b{x) + 6{x)^ =1 



x=0 



This leads to a quadratic equation for B 



T. 



[B'^f + 2UB^ + V" = , 



where 



N-l 



u 



V 



a:=0 

— (2 VpM + m - K^)) - K 



(74) 

(75) 

(76) 
(77) 



bmce ^/p{x) - h{x) > 0, using the inequalities (fTOj) and (f?^ together with the bound 
\/p{x) < 1/ y/rjN we obtain 



2" 3a/2 2"** 

-2T^= < [/ < ^=(l + 2T2"''/2), 



VV ~ ~ ^/T]N 

2~ii /2~"^ \ 2~3a/2 



-4T f h2""/2^ < r < 4T 

r]N\2T J - 



7]N 



■1 + r2-'^/2) _ 



(78) 
(79) 



As mentioned earlier, we assume that T is chosen to be the smallest integer for which 

<2-\ a 



2T2 



0) 



The above bounds can then be simplified as follows. 

2-a 



\U\ < 2 



\V\ < 8T. 



2-3a/2 

r]N 



(81) 
(82) 



11 



The value of therefore satisfies the bound 

< \u\ + ^/u^ + \v\ 



< 2 4= + x/9T2-3«/2/(r^iV) 

y/rjN V 

2-3a/4 

< aVt . (83) 
Using this bound together with (f7^ we obtain from Eq. (fTTjl the result 

I (^r I I > 1 - - (4yr2-3a/4 + + 2T2-''/2 j 

2-a/2 

> 1 - 3T . (84) 

V 

Directly from Eq. (|S^ we obtain the upper bound on the failure probability, 

= (2'^iV - N)\B^\' < 16T ^ . (85) 

3.3 Adding phases 

The state resulting from the measurement of the a auxiliary qubits has real ampli- 
tudes, i.e., it is of the form 

\^^) = "^^/p(x)\x) . (86) 

X 

The final stage of the algorithm, see Eq. (fT^ . turns into the final state |\E'), which 
can be written as 

1^) = X] V^(^) exp[2m(f){x)] \x) , (87) 

X 

where the 4>{x) are T'-bit approximations to the target phases 0(x), i.e., 

|0(x) -0(x)| < 2"^' . (88) 
Putting everything together, we find 

1(^1^)1 = I a/p(x)p(x) exp[27ri(0(x) - 0(x))] 

X 

> a/p(x)p(x) cos[0(a;) — (f>{x)] 

X 

> J]VR^)M(l-[0(a;)-0(x)]V2) 

X 

> J]V^0S^(i-2-'^'-^) 

= (1-2-2^'-^) 

> 1_3T^ (1-2-2^'-^), (89) 



V 

which is the required lower bound for the overall fidelity of the prepared state. 
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